— 


SCWOANDAURWNE 


. Introduction 
. Preview of Comm System Components 
. Preliminaries 


The Communications Channel 


. Analog Communication 

. Random Signals and Noise 

. Complex-Baseband Equivalent Channel 

. Digital Communication 

. Discrete-Time Implementation of Digital Communication 
. Error Analysis of Digital Communications 


Introduction 

This module presents introductory material on data communication, 
focusing on a high-level summary/comparison of analog and digital 
communication. 


Goal: Transmit a message from one location to another. 


When message is a... 

continuous waveform, use analog comm (e.g., FM radio), 

sequence of numbers, use digital comm (e.g., mp3 file), 

but keep in mind that the sequence of numbers might represent a sampling 
of a continuous waveform (as in the case of digital audio). 


Typical communication media: 


twisted pair wire (e.g., telephone,) 

coaxial cable (e.g., TV 4 p, datap) 

fiber optic cable (e.g., ethernet) 

EM waves (e.g., cellular phones 4p, WiFip, TV 4,p) 
water waves (e.g., underwater network, p) 


power lines 4p 
compact discp 
hard drivep 


magnetic tape 4 p 


where 4 = analog and p = digital. 


Note that, whether the message signal is discrete-time or continuous-time, 
the transmitted signal is continuous-time! 


Analog Communication 


analog transmitted received recovered 
message message signal de- message 
eS | modulator channel . le eee 


e Perfect recovery is impossible in the presence of noise! 


Digital Communication 


binary symbol baseband 
message | coder/ sequence pulse message 
———> J_____*____ 


---010110--| Mapper |) ...3-1,1,3... shaper Lua 


baseband passband received received 

message message passband de- baseband 
|__| modulator channel. _-————_» 

modulator 

received received equalized recovered 

baseband sequence . sequence de-mapper| message 

|| sampler equalizer --———_» ne 

|--3.1,0.9,1.1,-2.9-.| /decoder |. .oi0110... 


e A digital message is converted to an analog message coding and pulse- 
shaping, and then transmitted using analog modulation. To recover the 
message, the received signal is demodulated, sampled, and digitally 
processed. 

e Perfect recovery is possible even in the presence of noise! 


Preview of Comm System Components 
This module gives a brief description of important communication system 
components: the modulator, the coder/mapper, and the pulse shaper. 


The Modulator: 


e Translates “baseband” analog message to “passband”: 


where f, is the “carrier frequency.” 
e There are two principal motivations for doing this: 


1. Often we want to communicate several messages simultaneously 
(e.g., TV, radio, voice). It's difficult or impossible to do this if 
they overlap in frequency! 

2. Wireless EM transmission/reception is much easier at higher 
frequencies, since need antenna length > are A= - is 


wavelength and c=3e8 m/s speed of light.) 


system transmission band A/10 
VHF (TV) 30-300 MHz 1—0.1m 
UHF (TV) 0.3-3 GHz 10-1 cm 
cellular 824-960 MHz 3.cm 


WiFi 2.4 GHz 1cm 


Notice that practical antenna length determines where different 
signal types can be transmitted. 


The Coder/Mapper: 


e Coder transforms sequence of message bits into an error-resiliant 
sequence of coded bits. 

e Mapper transforms coded bits into discrete “symbols.” Ex: If the 
“symbol alphabet” is {—3, —1, 1, 3} and the symbol mapping is 


bits symbol 
00 3 

01 -1 

10 dL 

11 -3 


then ASCII text would be transmitted via 


letter ASCII code symbol sequence 
a 01 10 00 01 -1 i) -3 -1 


b 01 10 00 10 -1 1 -3 1 


Pulse Shaper: 


¢ Converts symbol sequence into a continuous waveform. 

e In linear modulation schemes, the time-n symbol a[n] scales a nT- 
delayed version of pulse g(t): 
Equation: 


m(t) = S\a[nlg(t — nT)” baseband message” 
T = ”symbol period” 


Ex: Say symbol sequence is [1, 3, —1, 1, 3]. Then 
, anjg(t — nT) for n =0,...,4 
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Preliminaries 

The module provides a review of the background concepts needed for the 
study of analog and digital communications: the Fourier transform, various 
definitions of bandwith, the Dirac delta, frequency-domain representation 
of sinusoids, frequency-domain plotting in Matlab, linear time invariant 
(LTT) systems, linear filtering, lowpass filters, and Matlab design of 
lowpass filters. 


Fourier Transform (FT) 


Definition: 
Equation: 


wi) = [ “ w(e Pd = F{w(t)} 


w(t) = f weap = F {WIN} 


Properties: 


e Linearity: F {cywy (t) + cow2 (t)} = c1Wi (f) + coW2 (f). 
conjugate symmetricW (f) 


e Real-valued w(t) => . 
ie |W (f)|symmetric aroundf = 0 


“Bandwidth” 


bandpass signal: 


5 a ef 


win) 
| a a, ae ae S-=== 0: dB 
ie ane Wate 
\ {1 4-power BW || 
| (= 
| 99%-power BW ! 
Ide = 


\ absolute BW \ . f 


Dirac Delta (or “continuous impulse”) 6(-) 


An infinitely tall and thin waveform with unit area: 


limit 
that's often used to “kick” a system and see how it responds 
Key properties: 


CO 


2. Time-domain impulse d(¢) has a flat spectrum: 
Equation: 


1. Sifting: if aia a. 


FAd(t)}= [- d(tje ?""dt = 1 (for all f). 


3. Freq-domain impulse 6(f) corresponds to a DC waveform: 
Equation: 


Ets ( f= a 6(fje?"Ftdf = 1 (for all t). 


Frequency-Domain Representation of Sinusoids 


Notice from the sifting property that 
Equation: 


ep / ” 8(f — folePVtdf = oP het, 


Thus, Euler's equations 


Equation: 
cos (afoot) = FePrfot + Le Pnfot 
sin (27 fot) = mer _ eee 


and the Fourier transform pair e??"/ot <> 6(f — f.) imply that 
Equation: 


F {cos (2a fot)} = 7O(f— fo) + 75(f + fo) 
F {sin (Qrfot)} = 376(f — fo) — 379 (F + fo). 


Often we draw this as 


|F {cos(27 fot) }| = |F {sin(27 fot) }| 


i 
2 
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Frequency Domain via MATLAB 


Fourier transform requires evaluation of an integral. What do we do if we 
can't define/solve the integral? 


1. Generate (rate-=-) sampled signal in MATLAB. 


2. Plot magnitude of Discrete Fourier Transform (DFT) using 
plottf.m (from course webpage). 


Square-wave example: 


f = 10; a 
t_max = 2; i * 
Ts = 1/1000; aay 
t = 0:Ts* t_max: a eS eS 2 
x = . 
Sign(cos(2"pi*f*t)); ss 
plottf(x,Ts); ‘& 


Noise-wave example 
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t_max = 1; 

1s = 171000° 

xX = 
randn(i,t_max/Ts); 
DLOttr( x, Ts) 


Notice that pLottf.m only plots frequencies f € | sp a). 


Linear Time-Invariant (LTI) Systems 


An LTI system can be described by either its “impulse response” A(t) or its 
“frequency response” H(f) = F{h(t)}. 


. casa impulse 5) | tte | h(t) — 
omain system response 


in inl flat 1 H(f) frequency 
omain spectrum system ‘response 
Input/output relationships: 


e Time-domain: Convolution with impulse response h(t) 


x(t) —-{ h(t) + u(t) y(t) = A(t) * a(t) = : " hé= Denar 


e Freq-domain: Multiplication with freq response H(f) 


X(f [HAE YS)  Y(f) = H(A)X(f) 


Linear Filtering 


Freq-domain illustration of LPF, BPF, and HPF: 


ee UT Og | 
eer ed, a, 


Lowpass Filters 


Ideal non-causal LPF (using sinc (a) := —— 


1 if\sB oF 
0 |f|>B 


h(t) = 2Bsinc(2Bt) 


Ideal LPF with group-delay tf: 


ei |f|/<B ; 
Binh = — > h(t) = 2Bsinc (2B(t—-t 
(f) fh > B (t) inc (2B(t — to)) 


|\H(f)| 
tho. 
-B 0 BS 


A causal linear-phase LPF with group-delay t,: 


to XT to 


symmetry around center yields linear phase 


but MATLAB can give better causal linear-phase LPFs... 


In MATLAB, generate 7~-sampled LPF impulse response via 


it hms (te. [Oxy hS.a |) 1G, G,e, Ol) is, 


where... 


Lf+1 = impulse response length 


“| | {0, fp}, {fs, 1} = normalized freq pairs 
a 
0 fp fs 1 


Tory 1G, G}, {0, 0} = corresp. magnitude pairs 


The commands firpm and fir2 have the same interface, but yield 
slightly different results (often worse for our apps). 


In MATLAB, perform filtering on 7 -sampled signal X via 


Ve= is tilter(i,., <). or Y= Fs -eonvi(n, x); 


tax = 3; Ts = 1/1006; 
xX = randn(1,t_max/Ts); 


h = firls(100,[0,0.2,0.4,1], 
[4,1,0,0])/1s; 

y = Ts*filter(h,1,x); 
subp.lot(3,1,1); 
DLOtLT (x, Ts,’ TF’ ); 
ylabel(’|X(f)|") 
subplot(3,1,2); 
plottr(h, Ts,7? Ff" j; 
ylabel(’|H(f)|") 
subplot(3,1,3); 
plOTTT(y, Ts,’ TZ 
ylabel("|Y(T)1") 


Note:The routines firls, firpm, fir2 generate causal linear-phase 
filters with group delay = = samples. Thus, the filtered output y will be 
delayed by = samples relative to x . 


The Communications Channel 
This module gives a brief introduction to the communications channel, 
which is modeled here as a linear filtering plus an additive noise. 


The effects of signal propagation are usually modeled as: 


the communications channel 


transmitted signal w(t) | received signal 
| ISI . r(t) , | RUDI 
\ | tf % | ‘| Sf 
| | ’ s(t) +> H(f) (+) 4» r(t) ; 
= i} | 
| 


where 


H(f): linear filtering due to multipath propagation 
w(t): additive noise/interference. 


Noise/interference sources include: 


e electronic circuitry (“thermal/shot noise” or “quantization noise”); 
usually broadband, 

e other comm systems (“multi-access interference” or “co-channel 
interference”); broadband or narrowband. 


SNR can be improved with appropriate filtering at receiver: 


a P 
¢ signal-to-noise ratio (SNR) = === S"sY — <s 
noise energy En 


e roughly speaking, average SNR can be improved by filtering out the 
frequencies dominated by noise: 


before filtering: 


after filtering: 


ra geet 4a 
ff 
uf 


Filtering due to Multipath Propagation 


The signal may propagate along paths with different lengths: 


a(t) F x(t) 


Since different path lengths imply different path delays: 
Equation: 


x(t) = cys (t — 71) + cos (t — 7) +--+ + ens (t — Ty), 


which can be written as x(t) = s(t)*h(t) for 
Equation: 


h(r) = c16 (t — 71) + €26 (7 — 72) +--+ + end (7 — Ty). 


The result is an |H(f)| that varies with f, implying frequency-dependent 
signal attenuation: 


t_max = 1; 

TS 1/1000; 

c1 1; taul = 0.02; 
C2 -0.5; tau2 = 
0.08; 

c3 = 0.2; tau3 = 


07.10; 


h= 
zeros(1,t_max/Ts); 
h(tau1/Ts) = c1i/Ts; 
h(tau2/Ts) = c2/Ts; 
h(tau3/Ts) = c3/Ts; 
plottf(h,Ts); 
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Analog Communication 

This module describes basic analog modulation techniques, including 
amplitude modulation (AM) with suppressed carrier, AM with a pilot tone or 
carrier tone, quadrature AM (QAM), vestigial sideband modulation (VSB), 
and frequency modulation (FM). Various demodulation techniques are also 
discussed, including envelope detection and the discriminator. Application 
examples include NTSC television and FM radio (both mono and stereo). 


1. Amplitude modulation (AM) 

2. Quadrature amplitude modulation (QAM) 
3. Vestigial sideband modulation (VSB) 

4, Frequency modulation (FM) 


AM with “suppressed carrier” 


AM of real-valued message m(t) (e.g., music) is 


ees, 
mi{t) aa s(t) s(t) = m(t) cos(27 fet), 

(~) cos(2z ft) fe = carrier freq. 

Nw 


Euler's cos (27 fct) = 4 |e?" + e-Prset| then implies 
Equation: 


si) = i m (t) cos (27f,t)e ?" dt 


= = / 7 m (the 2™F-fodtde + cs i * m (t)e 22Ft felt de 
2 picks Ca hee 
1 1 


aM Uf - fe) en 9 MF + fe): 


4 | M( 
1 NM(A) KIS(A)| 2B 


—fe fe 


Because m(t) € R, know |M(f)| symmetric around f = 0, implying the 
AM transmitted spectrum below f, is redundant! This motivates the QAM and 
VSB modulation schemes... 


With f. known, AM demodulation can be accomplished by: 


r(t) a v(t) u(t) = LPF{r(t) - 2cos(2z fet) }. 


(~) 2 cos(27 ft) 


For a trivial noiseless channel, we have r(t) = s(t), so that 
Equation: 


u(t) = LPF{s(t)-2 cos (27f.t)} 


= LPF m(t)-2cos? (27f.t) 
1+cos(2n-2f,t) 
LPF {m (t) + m(t) cos (27 - 2f-t)} 
m(t), 


assuming a LPF with passband cutoff B, > W Hz and stopband cutoff 
B, <2f.—W Hz: 


—~2f-e —2f-+W —f. —B, —W W B,B, fe 2f--W 2f- 


Note that we've assumed perfectly synchronized oscillators! 


When the receiver oscillator has {freq,phase} offset {, y}: 
Equation: 


u(t) = LPF m(t)cos (27f,.t) -2 cos (27 (fe + y)t + 9) 
cos(2ryt+y) + cos(27(2f.+y)t+y) 
m (t)cos (27yt + y) 


time-varying attenuation! 


Note:a freq offset of \ = Ze : 


v m/s between transmitter and receiver. 


Hz can occur when there is relative velocity of 


AM with “pilot tone” or “carrier tone” 


It's common to include a pilot/carrier tone with frequency] f-: 
Equation: 
s(t) = m/(t) cos (27f,t) + A cos (27f-t) 
pilot /carrier tone 
= |m(t) + A] cos (27f.t) 
S(f) = yIM(f-f.)+M (f+f.)+Ad (f- fe) + Ad (f+ fe) 


KIS(f)| 


Advantage: aids receiver with carrier synchronization. 
Disadvantage: consumes transmission power. 


While modern systems choose A <max |m(t)|, many older systems use 

A >max |m(t)|, known as “large carrier AM,” allowing reception based on 
envelope detection: 

Equation: 


| 


v(t) = [LPF {|r(@|}-A 


2 


m(t) (with a trivial channel) 


where |-| can be easily implemented using a diode. 


The gain + above makes up for the loss incurred when LPFing the rectified 
signal: 


desired output level —— 


|- 


Jo* cos(2rfet)dt 2 
he 


LPF output level —— 


— 


< MATLAB code here >> 


Quadrature Amplitude Modulation (QAM) 


QAM is motivated by unwanted redundancy in the AM spectrum, which was 
symmetric around fc. 


QAM sends two real-valued signals {m, (t), Mq (t)} simultaneously, 
resulting in a non-symmetric spectrum. 


my(t) —» “in phase” 
cos(27 f.t) 


~ (> s(t) = m,(t) cos(27f,t) — ma(t) sin(27 ft) 


sin(27 f,t) 
ma (t) Smead “quadrature” 


QAM demodulation is accomplished by: 


2 cos( ca ft) 


2 oor fet) 


| LPF 


where the LPF specs are the same as in AM, i.e., passband edge B, > W Hz 
and stopband edge B, < 2f, — W Hz. 


For a trivial channel, we have r(t) = s(t), so that 
Equation: 


uz(t) = LPF{r(t)-2 cos (27f,t)} 


= LPF my,(t)2 cos” (27f,t) 


1+cos(47 f,t) 


—mMg (t)2 sin (27 f-t) cos (27 ft) 
sin(47 ft) 
= mr(t) 
vg(t) = LPF{-r (t)-2 sin (27f.t)} 


= LPF -—m,(t)2 cos (27f,t) sin (27f,t) 


sin(47 ft) 


+mg (t)2 sin? (27 ft) 
1—cos(47 f,t) 


= mg (t), 


assuming synchronized oscillators. 


When the oscillators are not synchronized, one gets coupling between the 
I&Q components as well as attenuation of each. Writing the I&Q signals in 
the “complex-baseband” form 

Equation: 


m(t) = mr(t) + jmg(t) 
w(t) = v(t) + jog () 


yields a much simpler description of QAM: 


QAM modulation QAM demodulation 


in(t) —G-| Re fo s(t) r(t) a ag lcd a(t) 


(~) el 2n fet (~) Je JI 2a fet 


4 |M(S)| 


A IRCA) 


A IMC - f)I 


A IS(P)I 


A IV(S)I 


Note:Re {u (t)} = + lu (t)+ur (t)| Sj + (U(f) ey (—f)]. 


We now verify the complex-baseband model for modulation: 
Equation: 


Re {im (t)e?7/t} 


Re {(mz (t) + ma (t))(cos (27 f-t) + j sin (2a f-t))} 
my (t) cos (27f-t) — mq (t) sin (27 ft) = s(t), 


as well as for demodulation (assuming r(t) = s(t)): 
Equation: 


w(t) = LPF{s(t)-2e 27} 
LPF { (my (t) cos (2afet) — mg (t) sin (2rfet)) -26e Pe} 
— LPF{m, (t) (effet +e Frhet) o— Pact 
—mMg (t) (Je pam fet _ jePrtet) e Pafty 
= LPF {mar (t)(1 +e) — ma GEM — 3)} 
= myz(t) + jm (t). 


The convenience of complex-baseband results in widespread use of complex- 
valued signals for comm systems! 


Note:To get the complex baseband formulation for AM, we simply set 
Mg (t) = 0 and m;(t) = m(t). 


Vestigial Sideband Modulation (VSB) 

VSB is another way to restore regain the spectral efficiency lost in AM. It's 
used to transmit North American terrestrial TV, both analog (NTSC) and 
digital (ATSC) formats. 


Like AM, it can operate with or without a carrier tone. 


Basically, VSB suppresses most of the redundant AM spectrum by filtering it: 


i, band 
m(t) +(x -+| ma s(t) 


{_™N 
~~) cos(27 ft) 


{_™ . 
i) cos(27 ft) 


The passband VSB filter is a BPF C(f) where 
Equation: 


C(f—f.)+C(f+fe-)=2 for |f| < W, 


which implies its inside rolloff is symmetric around f = f;: 


~2f, —fe fe *fe 
IC(f + fe) +C(F — fe)| 


For VSB modulation, we have 
Equation: 


s(t) = m(t) cos (27f,t)*c (t) 
Sf) = SIM(f+ Ff.) +M(F- FIC UH). 


It turns out that VSB demod is identical to AM demod: 
Equation: 


u(t) = LPF{r (t)-2 cos (27f.t)} 
= LPF{s(t)-2 cos (27f.t)} (trivial channel) 


V(f) = LPF{S(f— fe) + S(f + fe)} 
1 


= FUPF{ IM (f) + MF 2f1C (Ff) 
HM (f+ 2fc) +M (AIC (F + fe) 

= M(NSIC(F-f) + OUR +f) 

=1 for fe[-W,W] 


= M(f). 


We note that the property 
Equation: 


FP {cos (2nf.t)e(t)} = F[C(F- fe) + OUF + fe) 


may be convenient, e.g., for testing whether a given filter c(t) satisfies the 
passband VSB criterion. VSB filtering can also be implemented at baseband 
using a complex-valued filter response @(t) which satisfies 

Equation: 


C(f)+C" (-f)=2 for |f| < W, 
generating the complex-baseband message signal ™ (t). The message can be 


recovered by simplying ignoring the imaginary part of the complex-baseband 
output 0 (t). 


VSB modulation VSB demodulation 


=o —s 0) 0) -@~UPE [Rel 


( ~) ) eI 2m fet (~ ©) Qe-I2*fet 
| IM(p) \ IRS) 
Ka \ = hs [ - Pi 


A IM(/I 


fh IR + fo) 


h IM(S — fe) 


=— f 
\ IS(P)| | IV(F)| 
f — — f 


Motivation: filtering at baseband is usually much cheaper than filtering at 
passband. 


Frequency Modulation (FM) 


While AM modulated the carrier amplitude, FM modulates the carrier 
frequency. 


t_max = 2.0; W=1; % message 
params hs 

Ts = 1/1000; t = 0:Ts:t_max; it 

m = sin(2*pi*w*t); m Message = = =.= 
signal 

fc = 20; % carrier 
freq 


D = 15; % FM mod 


index 

kf = D*W/max(abs(m)); % freq 
sensitivity 

S_am = m.*cos(2*pi*fc*t); 

s_fm = 
cos(2*pi*fc*t+2*pi*kf*cumsum(m)*Ts); 
subplot(3,1,1) 

plot(t,m), 

grid on; title('message'); 
subplot(3,1,2) 

plot(t,s_am); 

grid on; title('AM'); 
subplot(3,1,3) 

pLot(t;s-fm); 

grid on; title('FM'); 


In particular, FM modulates the real-valued message m(t) via 
Equation: 


s(t) =cos 2nf.t+ Qrke fim (r)dr 


o(t)” instantaneous modulation phase” 


where ke is called the “frequency-sensitivity factor.” Since the instantaneous 
modulation frequency 
Equation: 


ao) — 2k m (t) 


is a scaled version of the message m(t), it is fitting to call this scheme 
“frequency modulation.” 


Using the peak frequency deviation A; = ks max |m (t)|, the “modulation 
index” D is defined as 
Equation: 

As 

W < denominator is the single-sided BW of m/(t). 


Increasing D decreases spectral efficiency but increases robustness to 
noise/interference. 
Equation: 


D<«1 =: **‘narrowband FM”, 
D> 1 : ~**wideband FM’’. 


Carson's Rule approximates the FM passband signal-BW as 
Equation: 


BWog © 2(Ap+W) =2(D+1)W. 


Example: Mono FM radio: 


e Message signal filtered to freq interval [30,15k] Hz. 
e FCC limits A; < 75 kHz (channels 200 kHz apart). =D = e = 


FM stereo uses smaller D due to message spectrum: 


other optional services 
he n L+R 
sum =; (e.g., text, muzak) 


AM-modulated “diff” ut ———eeN 


} 
} 
j 
/ 
} 
| 


i 


kHz 


15 19 38 57 76 95 


There are various FM demodulators, but the “discriminator” is one of the best 
known. Recalling that 


Equation: 
dp(t 
© cos (6(¢)) =- 9 sin (9(t), 
we see that 
Equation: 
d 
ee (t) = cos (2nfet + Qk fim (r)dr) 


= —[2rf. + 2aksm (t)] sin (2nfet + Ir to m (r)dr) 


is a form of large-carrier AM (assuming f, > kp m (t)), which can be 
demodulated using an envelope detector as follows: 


d envelope 


< MATLAB code here > 


Random Signals and Noise 

This module describes the basics of random signals and noise as needed for an 
introductory course on data communications. First it introduces the notion of 
power spectral density (PSD) and defines "white" noise as a noise with a flat 
PSD. Next it describes the effect that linear filtering has on a random signal, 
in both time and frequency domains. For this, the autocorrelation function is 
introduced, and expectations are used to derive the fundamental properties. 


When designing a comm system, it is impossible to know exactly what the 
signal and noise waveforms will be. But usually we know average 
characteristics such as energy distribution across frequency. It is exactly these 
Statistics that are most often used for comm system design. 


We will consider random waveforms known as “zero-mean wide-sense 
stationary random processes.” Such an z(t) is completely completely 
characterized by its power spectral density S, (f), or average signal power 
versus frequency f. The technical definition of PSD accounts for the fact that 
power is defined as energy per unit time: 

Equation: 


2 


1 : . 
as j2nft 
Sx (Ff) = im oT E [avec dt 


Above, E {-} denotes expectation, i.e., statistical average. 


A common random waveform is “white noise.” Saying that w(t) is white is 
equivalent to saying that its PSD is flat: 
Equation: 


Sw (f) = No for all f. 


Broadband noise is often modelled as white random noise. A fundamentally 
important question is: What does filtering do to the power spectrum of a 
signal? The answer comes with the aid of the autocorrelation function R, (rT) 
, defined as 

Equation: 


R, (7):=E {2(t)x" (¢-7)} 


and having the important property 
Equation: 


R, (7) © Sz (f). 


Note that, for white noise w(t) with PSD N,, we have 
Equation: 


S,(f)=N, => R, (tT) =N,6(1) 
— £E {w (t)w" (t — r)} = N,6(r). 


Here we see that w(t;) and w (t2)|t2 A t1 are "uncorrelated" because 

E {w (t;)w- (t2) } = 0, and that w(t) has average energy 

Ey = E {|w (t)|"} = 00, which may be surprising. We could have found the 
same via Sy = f°. Sw (f)df. 


Now say that white noise w(t) is filtered with non-random A(t), yielding 
output y(t) = f°, w(q)h (t — q)dg. We can find Sy (f) by first finding 
R, (7) and then taking the FT. Recall: 

Equation: 


Ry (rt) =E {y(t)y" (t= 7)}. 


Plugging in the expression for y(t), we find 
Equation: 


me) 
Ned 
a 
A 

II 


E if (qi)h (t — qi)dqy [ w (q2)h° CTS an)daa 


—0oO 


B{f fw (q1)w" (q2)h(t—qi)h° (t—7 - wn)danda 


Because averaging is a linear operation, and because the h(-) terms are fixed 
(non-random) quantities, we can write 
Equation: 


= [fb {wlayw" ed} me ayn 7 ~ an)darder 


Since w(t) is white, 
Equation: 


E {w (a1)e0" (a2) } =E {w (an)e" (a — (a1 ~ a2) } = Nod (ar ~ 92), 


which allows use of the sifting property on the inner integral: 
Equation: 


Ry, (7) 


f a N05 (q1 — g2)h(t —q)h (t — 7 — q2)daidae 


No te! (t — qo)h (t — 7 — qa)dqe 


h(q)h’ (q—7)(—dq) using g:=t- @ 


| 
Z 


h(qh (q—7)dq. 


| 
2 
ime 


To summarize, the autocorrelation of filtered white noise is 
Equation: 


R, (7) = No [x (q)h’ (q—T)dq. 


Having R, (7), the final step is finding S,(f) = A {R, (7)}: 
Equation: 


SG = [ . R, (r)e Pdr 


[| i0 [mean (@— rag] eter 

= N, [on (qe Ia ie h* (q—- nett Nar| dq 
= N, [on (qje P?7F9 if h” (hear dq 

= No fon (qje "dq if h (je*hat ) 


= NH (f)H*(f) = NolH (A)? 


To summarize, the power spectrum of filtered white noise is 
Equation: 


Sy (f) = Nol (f)/?. 


More generally, the power spectrum of a filtered random process x(t) can be 
shown to be 
Equation: 


Sy (f) = Sz (f)|H (AY, 


which is quite intuitive. 


Complex-Baseband Equivalent Channel 
In this module, the "complex-baseband equivalent channel" is derived from 
the three-part cascade of modulation, wideband channel, and demodulation. 


Linear communication schemes (e.g., AM, QAM, VSB) can all be 
represented (using complex-baseband mod/demod) as: 


modulation a ; demodulation 
_ a s(t) DI ; 
nt) mL ne opm fr 8 
~ | ¥ < 
( (~) el2T fet | w(t) ! (~ ) De—I2a fet 


It turns out that this diagram can be greatly simplified... 


First, consider the signal path on its own: 


s(t) 
so p-ale le e 


~) ef 2m Set ( (~) ) QeI2R fet 


Since s(t) is a bandpass signal, we can replace the wideband channel 
response h(t) with its bandpass equivalent hpp (t): 


KIS, ISCO 


...because filtering 
s(t) with h(t) is 
equivalent to 
filtering s(t) 

with h,,(t): 


j |e) 


- f 
\ |S(f) Hoo( F)| 


Then, notice that 
Equation: 
[s (t)*hpp (t)|2e 2274 8 (T)hpp (t — T)dr - 2e P74 


s (7) 26° Phy, (t—r)e Pedr 


| 


| 


(26 2Ne|* | hop (jen) ; 


which means we can rewrite the block diagram as 


- s(t) - : 7 
4 ei 27 fet 4 Je —I27 fet 


We can now reverse the order of the LPF and hp, (t)e 7?""“ (since both are 


LTI systems), giving 
hy, (t)eJ27 F<! me H(t) 


Since mod/demod is transparent (with synched oscillators), it can be 
removed, simplifying the block diagram to 


Now, since m (t) is bandlimited to W Hz, there is no need to model the left 
component of Hip (f + fe) = F {hpp (the PM}: 


j |H(£)| 


‘pA f + fe)| {AA 
> f ae i 


—fe fe ake. 0 —B 0B 


Replacing hy, (t)e 74 with the complex-baseband response h(t) gives 
the “complex-baseband equivalent” signal path: 


m(t) Us(t) 


The spectrums above show that hp, (t) =Re {h (t) - deirrit\, 
Next consider the noise path on it's own: 
w(t) (| LPF = on (t) 


aap, 2 
Saal %e j2n fet 


iW (f)| 
Wy y y (Wry, 
ar aa J 


\IW(f + fo)| 
PARI TAM, 


From the diagram, @,, (t) is a baseband version of the bandpass noise 
spectrum that occupies the frequency range f € [f. — Bs, fe + Bs]. 


Since @, (t) is complex-valued, V;, (f) is non-symmetric. 


Say that w(t) is real-valued white noise with power spectral density (PSD) 
Sw (f) = No. Since Sy (f) is constant over all f, the PSD of the complex 
noise @, (t) will be constant over the LPF passband, i.e., f € [—Bp, B,|: 


Si, (Ff) 


-B, -B, _W W B, B, f 


A well-designed communications receiver will suppress all energy outside 
the signal bandwidth W, since it is purely noise. Given that the noise 
spectrum outside f € [—-W, W] will get totally suppressed, it doesn't 
matter how we model it! Thus, we choose to replace the lowpass complex 


noise U,, (t) with something simpler to describe: white complex noise w (t) 
with PSD S.,(f) = No: 


-B, -B, _W WB, B, J 


We'll refer to @ (t) as "complex baseband equivalent" noise. 


Putting the signal and noise paths together, we arrive at the complex 
baseband equivalent channel model: 


' complex baseband 


| 
1 equivalent channel 


! Ny, (t) = Re{h(t) - 2672" et 
ae foe ayy, nf = Rati) -20mr 


w(t) = complex-valued white noise 


complex baseband bandpass equiv complex baseband 
modulation channel demodulation 


m(t) (x) 


1/Q modulation i : Q demodulation 


m,(t) — x Our = u(t 


cas(27 f.t) 
p ms" 
vy, 


: = \ 
sin(2r f.t) as - 
Ma(t) —+ —»(x) (t x) VQ(t) 


The diagrams above should convince you of the utility of the complex- 
baseband representation in simplifying the system model! 


Digital Communication 

This module outlines the basics of digital communication through a white- 
noise channel. Transmission consists of pulse shaping and modulation, 
while reception consists of demodulation, filtering, and sampling. The 
combined pulse-shaping and filtering designs which prevent inter-symbol 
interference (ISI) are identified as those satisfying the Nyquist criterion, 
which is examined in both the time and frequency domains, and the raised- 
cosine (combined) pulse is given as a common example. The Cauchy- 
Schwarz inequality is then used to show the matched filtering maximizes 
the signal-to-noise ratio. Finally, the square-root raised cosine pulse/filter is 
given as an example of matched filtering that satisfies the Nyquist criterion. 


Transmission consists of 


1. pulse shaping: m (t) = So, a[n]g(t — nT), 


Reception consists of 


1. demodulation:  (¢) = LPF {2r (t)e #277} 
2. filtering: y(t) = 6 (t)*q (t), 
3. sampling: y[m] = y(mT). 


digital modulation digital demodulation 
m(t) s(t) \ 7 r(t) ne 
fof}! | Yolo 
& 
u 
. nT ~ el2t fet ~ %e jn fet 
analog QAM mod analog QAM demod 


Building on analog QAM mod/demod components, digital mod adds pulse 
shaping & demod adds filtering/sampling. 


Simplifying via the complex-baseband equivalent channel: 


a le Oe 


a) nT Ww 


Transmitter pulse shaping is used to convert the symbol sequence {a|n]} 
into the continuous message ™m (t): 
Equation: 


in (t) 


S- aln|jg(t—nT) baseband message” 


T = ”symbol period” 


Thus, 7 (t) can be seen to be a superposition of scaled and time-shifted 
copies of the pulse waveform g(t). 


Example, if the symbol sequence [a[0], a[1], a[2], a[3], a[4]] equals 
[1, 3, —1, 1, 3], then the square pulse g(t) shown below left yields the 
message m (t) shown below right. 


, a[n|g(t — nT) forn = 0,...,4 


T 2T 37 47 5T 6T 


g(t —T): + 
~ t ~ t 
g(t — 2T): 4 
>t 
g(t —3T): + 
oa * 
g(t —4T): 4 
~ t 


Receiver filtering (via g(t)) has two goals: 


1. noise suppression (i.e., SNR improvement), 
2. inter-symbol interference (ISI) prevention. 


Noise suppression was briefly discussed in Preliminaries and will soon be 
revisited in more detail. Next we describe ISI. 


Realize that, in the ideal digital comm system, the n*” output y[n] would 
simply equal the n“” input a[n]. But in practice, y[n] can be corrupted by 


interference from the other symbols {a [m]},,.2,» known as “inter-symbol 


interference,” and noise. 


ISI-prevention for the noiseless trivial channel 


Consider the idealized system 


m(t 
a|n| s(t) ©) 


@ nT 


Equation: 


sey. = / a(r)i(t—r)dr for m(t) => ang (t— nT) 


q(r)g(t — nT — 7r)dr 
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g 
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| 
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Thus, the idealized system can be re-drawn as 


where 
Equation: 


y[m] = y(mT) = } Ja [nlp(mT — nT) = Sa [n]p((m — n)T). 


To make y|m] = alm] (i.e., prevent ISI), we need 


1 m=0 
0 m0 


p(mT) = 
—3T-2T -—T 0 T 2T 37 
which is known as the “Nyquist Criterion.” This criterion can be simply 


stated as p(mT’) = d|m] using 
Equation: 


a|n] = », a|m]é[n — ml] ’sifting property.” 


mMm=— Oo 


Examples of Nyquist, and non-Nyquist, combined-pulses p(t) for 
[a [0], a[1], a[2], a[3), a[4]] = [1, 3, -1, 1,3): 


, ain|p(t — nT) for n = 0,...,4 


t 
t t 
t 
3 ye e 
t 
1 ° 
x z aT 47 5ST tT 


rawvanre a{n|p(t — nT) for n = 0,...,4 


p(t —2T)- 


n|p(t — nT) 


p(t —3T):- AN 


y(t) =>, 


T 27 3 4T ST 6T 
e 


There is an interesting frequency-domain interpretation. Since 
Equation: 


= 3 (1-5) 4 S 5(t—mT), 


k=—00 mMm=— OOo 


we can see that 
Equation: 


Pte > o(f-F) & pe: Yo see—mry. 
k=—0o mMm=—oo 
1 k = 
Z > Ps 7) DY pla (emery 


So, the time-domain Nyquist criterion p(mT') = 6[{m] implies 
Equation: 


= 3 P(f- 5) & te), 


k=—00 


which in turn implies 


In other words, the superposition of { me ( f- #) \ pez, Must sum to one. 
This frequency-domain version of the Nyquist Criterion will soon come in 
handy... 


A popular choice of combined pulse p(t) = g(t)*q(t) is the “raised-cosine 
pulse” with rolloff parameter a € [0, 1}: 
Equation: 


cos (zat /T) sin (7x 


prc(t) = 1 (at/T) sinc (t/T), sinc (x) := = 
fi i. 

Pao(fy = « Tees? (2 (7/2). << fie. 
0 Cre) < |f| 


time (T=1) 


-1 -0.8 -0.6 -04 -0.2 0 0.2 04 0.6 0.8 1 
frequency (T=1) 


So, we now know how to design the combined pulse p(t). 


But what about the individual pulses g(t) and q(t)? 


Maximizing SNR for ISI-free Pulses in White Noise 


Now let's bring the noise back into consideration. Given 


alr] —=| a(t) FR 


Equation: 


E {a (t)® (t— r)} = Nod(r) (complex white noise) 


we want a {g(t), g(t)} pair that maximizes the SNR of y|m]. 


Separating the noise and signal contributions to y|m] via 


w(t) = 


an] 50) 


p(t) 


the SNR can be written 
Equation: 


& _ E {lys(mll?} 
6, E {lyn [m[P} 


where E, and E,, are average signal and noise energies. 


Here, we treat both @# (¢) and a/n] as random, implying that y, [7m] and 
Yn |m] are both random. 


Notice that, with an ISI-free combined pulse p(t), we get 
Equation: 


ys(m] = ) a[n|p((m—n)T) = a[m)p(0) 


3 

—~ 
=) 

— 
| 


[ a@ao-nar 


(@) 


so that 
Equation: 


6 = Eflyslml|?} =E{{alm| [ a(r)g(—rar 


Oo 


J 


) 


— | aea(-nar 


where o,7 denotes average symbol energy. Next, notice that 
Equation: 


ia = aecany / Wa Grease 


Oo 


so that 


Equation: 


én 


E {yn [m|?} =e {| [90 B ont ~ tf 


= E if q(t) (mT — nar q (r')® (mT - r')dr'} 


= [ q(T) [- q (r'\E {a (mT — 1)@ (mT — 7’) \ar'dr 
—0oo —0oo . a es 
Noé (r' — 7) 


= / la(7) 2dr. 


Putting these together, we find 
Equation: 


o? | fr 7)g (—r)dr|" | 


SNR = 
7 i ze ) [Pdr 


Cauchy-Schwarz says 
Equation: 


ee (r)dr|” < f%lb(r)Pdr- f%. le (7) Par 


with equality iff b(7) = Ke (r) for any K, 


which implies 
Equation: 


SNR < a [i (—r)|?d 
—T 1 
NG it gare 


with equality iff g(r) = Kg (—r) for any K. 


Noting that SNR doesn't depend on K, we choose K = 1. Thus, given pulse 
g(t), the SNR-maximizing receiver filter is 
Equation: 


* 


q(t) =g (—T) knownasa”matched filter”. 


We can write this in the frequency domain as 


Equation: 
a) = fate tran = fo gf (re Ptar 
= [os (t)e?"dt = ic jetta) 
ae) : 


Summary: For SNR-maximizing ISI-free pulses, we need 


1. G(f)Q(F 


) = P(f) satisfying the Nyquist criterion, 
2.Q(f) =G 


); 


which together imply |G (f) |? must satisfy the Nyquist criterion. 


One option is G (f) = \/ Pra (f), since Pac (f) was Nyquist. We call this 
the “square-root raised cosine” (SRRC) pulse. Working out the details of 
aes {Ggrrc (f)}, we find 

Equation: 


(1—a) sinc (4(1—a))  4acos (r£(1+a)) 


1 — (4a4)* n(1 — (404)”) 


gsrerc (t) = 


—1 -0.8 -0.6 -0.4 -0.2 0 0.2 04 0.6 0.8 1 
frequency (T=1) 


At the receiver, we would use gsrrc (t) = ggrrc (—t) = gsrrc (t); the 
latter equality is due to gsrrc (t) being real and symmetric. 


Discrete-Time Implementation of Digital Communication 

In this module, we derive the discrete-time implementation of a (baseband- 
equivalent) digital communication system, which is useful for Matlab 
simulation and error-rate analysis. To do this, we use the fundamental DSP 
tool known as "sinc reconstruction", which says that a bandlimited 
waveform can be represented exactly by a sum of shifted sinc pulses. Sinc 
reconstruction is separately applied to the transmission pulse, the reception 
pulse, and the baseband equivalent channel impulse response, after which 
those three steps can be represented by discrete-time filtering operations. 
For this, we show that it is essential to sample the filters at least twice as 
fast as the symbol rate. 


Digital implementation of transmitter pulse-shaping and receiver filtering is 
much more practical than analog. First, recall the analog implementation of 
our digital communication system: from the module Digital 
Communication. 


a{n] — s(t) — in KO 
bi 


(Ly) nT Ww t) 


To proceed further, we need an important DSP concept called “sinc 
reconstruction”: 
Equation: 


1 
If waveform z(t) is bandlimited to op Hz: then 


Equation: 


ne 3 zx |n] sinc (= -nt.)) for x(n] = x (nT;). 


nN=—0o 2 


In other words, a bandlimited waveform can be reconstructed from its 
samples via sinc pulse shaping. 


x(t) = > 5, 2[n] sinc(=(t —nT’)) 


Discrete-Time Pulse-Shaping 


Applying +, -sampling and reconstruction to g(t) (where the SRRC pulse 
bandwidth a requires the use of P > 2), 


Equation: 
g(7T) = — sinc (4 (7 — 15)) for gl] =g (1) 
m = n\g (t — nT) 


>, 9 [I] sinc (4 (t — nT -14)) 

9 [k— nP] sinc (E(t -k4)) viak=nP-I 
x, >>, a [n|g [k — nP] sinc (F(t = k<)) 

=m |k] 


a 


The sequence m [k], a weighted sum of P-shifted pulses g[k], can be 
generated by P-upsampling a[n] (i.e., inserting P—1 zeros between every 


pair of samples) and filtering with g|k}: 


al2] 


me alk AY sine( Be) 


Note:sinc-pulse shaping = digital-to-analog conversion (DAC). 


Discrete-Time Receiver Filtering 


Applying +, -sampling and reconstruction to bandlimited g(r): 
Equation: 


a7) = LE all] sine (= (7 —I5)) for a[l] =a (Us) 


where again we need P > 2, yields 
Equation: 


yk] = (hE) = SX a(7)6 (kB -7)dr 
= [° LE wall sine (E(r-1E)) o(6 —7)dr 
-> a[ll f°, sinc wie 5((k B= "ar 
pee t)*O()},_@ ye =Olk-1 
- > aloe — glk*o ry, 


from which y|m] is obtained by keeping only every P* sample: 


{sinc( rt) } 


i.e., downsampling. Here sinc (+ t) does anti-alias filtering. 


Discrete-Time Complex-Baseband Channel 


Finally, we derive a discrete-time representation of the channel between 
m [k] and 0 |k}: 


Using @ |k] to refer to the noise component of @ [kK], it can be seen from the 
block diagram that 
Equation: 


wk] = f° @(r) sine (£(k45 — 1)) dr 


To model the signal component of @ [k], realize that ™ |k] is effectively 
pulse-shaped by sinc (+t) *h (t)*sine (+t). But since the frequency 


response of sinc (+t) has a flat gain of 4 over the signal bandwidth, and 


thus the bandwidth of h (t), 
Equation: 


sinc (4t) *h, (t)*sinc (+t) es 


| 
Ps 
— 
oH 
— 


2 ~ 
So, with £ -sampling and reconstruction of ir h(t), i-e., 
Equation: 


EeR() = OE hldsine (F(t-éh)) for ble] = Feb (EF) 


1=—0O P 


we can write 0 [k] as 
Equation: 


o[k] = Gk] +d, Ml] Fh (kG - 15) 
= O[k] +32, m[l] 4, h[dsine (k — 1-1) 
d[k —l—i] 


yielding the discrete-time channel 


w{k] 


nk] | hk] Rete of] 


Merging the discrete-time channel with the discrete-time modulator and 
demodulator yields 


wk] 
] 


known as the “fractionally sampled” system model. This model is very 
convenient for MATLAB simulation and acts as a foundation for further 
analysis. 


Error Analysis of Digital Communications 

In this module, we first introduce the eye diagram and constellation diagram 
as qualitative ways of evaluating the symbol error probability of a digital 
communication system. We discuss various symbol alphabets, such as QAM, 
PAM, and PSK, and their associated decision regions. Finally, we derive the 
symbol error probability for PAM and QAM in additive white Gaussian 
noise, using the Q and erfc functions, and discuss Gray coding. 


Recall the figure below, from the module Discrete-Time Implementation of 
Digital Communication, When the channel is trivial and noiseless and the 
pulses satisfy the Nyquist criterion (i.e., g|/k]*q|k] = 6[k]), the digital comm 
system will work perfectly, yielding y|n] = a[n]. 

wk] 


| 


In practice, however, 


e the pulses g[k] and q|k] will be truncated to finite length, 


e the channel will not be trivial (i.e., h [k] 4 6 [k]), and 
¢ the channel will not be noiseless (i.e., w [k] A 0), 


leading to y[n| 4 a[n], in which case we must infer the value of a[n] from 


the received samples {y {m]}°°__.. For now, we consider using only the 


single sample y|n] to infer a[n]. 
Key question: What are the mechanisms by which errors are made? 


To better understand error behavior, we can plot the “eye diagram” or the 
“constellation diagram” and calculate the symbol error rate (SER). 


Eye Diagrams 


Usually used when a|n] € R, the eye diagram is a plot which superimposes 
T-second segments of Re {y(t)} over the time intervals 


t € [nT — 4,nT + +) for many integers n. 
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diagram 


bebe 


Re{y(t)} , 
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In MATLAB, the eye diagram can be made by superimposing P-sample 
segments of Re {y+ [k]} corresponding to intervals 
ke {nP ~ £, ey nP+ fo} for many n. (Usually P > 8.) 


As the eye “closes,” 
decisions get more 
unreliable: 


When the eye is 
“open,” decisions 
will be reliable: 


If aln] € C, eye diagrams can be plotted for both the “I” and “Q” channels 
using Re {y(t)} and Im {y(t)}, respectively. MATLAB for digital 
mod/demod with eye diagram: 


"% design SRRC ae 
P = 16; % oversampling factor ¥4¥.. 
alpha = 0.5; % SRRC rolloff i 
param : 
D = 2; % truncation to [-DT,DT] + 
g = srrc(D,alpha,P); % SRRC ast)! 


SE) ft hawt 


pulse z| ) 
Ng = length(g); aioe By La yg 
% generate symbols a 
N = 100; % # symbols (ncn at 


M= 4; % alphabet size 

sig2a = 1; % symbol variance 
a = pam(N,M,sSig2a); % symbol 
sequence 

% pulse-amplitude modulate 
a_up = zeros(1i,N*P); . 
a_up(1:P:end) = a; % upsampled | 
symbols 

m = conv(a_up,g); % PAM 

% matched-filter demodulate 
y_up = conv(m,g); % use SRRC 
again 

% remove causal filtering delay 
left = -P/2; % left edge of eye 
diagram 

k = left+[1:P*N]; % desired 
time indices 

dly = (Ng-1)/2+(Ng-1)/2;% delay 
due to pulses 

y_up = y_up(k+dly); % remove 


delay 

y = y_up([1-left:P:P*N]);% 
downsample 

% plot received signal 
figure(1) 


plot(left/P+[0:P*N- 
1 J7P; Yup, Fr", (O2Na1 Ve 3); 
axis( ‘tight’ ) 


title([’SRRC 
(\alpha=‘,num2str(alpha),... 
‘) truncated to 

Xp’ .num2str(D)..7T? |.) 
Xlabel([ ‘symbol index’ ]) 

% plot eye diagram 

figure(2) 

Y_up = reshape(y_up,P,N); % 
extract segments 
plot(left/P+[0:P-1]/P,Y_up) % 
Superimpose them 
title([’SRRC 
(\alpha=‘,num2str(alpha),... 
‘) truncated to 

\pm’, num2str(D),’T’]) 
Xlabel([‘relative symbol 
index’ ]) 


Constellation Diagrams 


The constellation diagram is a plot of Im {y[n]} vs. Re {y[n]} for many 
integers n. When the comm system is working well, the points cluster around 
the symbol alphabet values: 


05 


=“ “05 0 os 1 15 


Recall that y[n] € C due to the complex-baseband channel model, regardless 
of whether a[n] € Ror an] € C. 


Sometimes it is instructive to superimpose a plot of Im {y; |[n]} vs. 
Re {y+ |n]}, which approximates the trajectory of y(t) in the complex plane: 


MATLAB for digital mod/demod with constellation diagram: 


% design SRRC 
P = 8; % oversampling factor 
alpha = 0.5; % SRRC rolloff param 
D = 2; % truncation to [-DT,DT] 
g = srrc(D,alpha,P); % SRRC pulse 
Ng = length(g),; 
% generate symbols 
N = 100; % # symbols 
M= 4; % (sqrt) alphabet size 
sig2a = 1; % symbol variance 
a = gam(N,M,sig2a); % symbol sequence 
% pulse-amplitude modulate 
a_up = zeros(1i,N*P); 
a_up(1:P:end) = a; % upsampled symbols 
m = conv(a_up,g); % PAM 
% matched-filter demodulate 
y_up = conv(m,g); % use SRRC again 
% remove causal filtering delay 
= [1:P*N]; % desired time indices 
dly = (Ng-1)/2+(Ng-1)/2;% delay due to pulses 
y_up = y_up(k+dly); % remove delay 
y = y_up(1:P:end); % downsample 
% plot received signal 
figure(1) 
ae P*N-1]/P, real(y_up), imag (y_ UD) ners 
r’, [O:N-1],real(y),imag(y),*."); 
Xlabel( ‘symbol index’); 
ylabel(‘I’); zlabel(‘Q’); 
view(20,30); % to see full trajectory 
%also try view(0,90), view(0,0), view(90,0) 
% plot constellation diagram 
figure(2) 
plot(real(y_up),imag(y_up),’y’,real(y),imag(y),’.’); 
Xlabel(‘I’); ylabel(‘Q’); 
title([‘SRRC (\alpha=‘,num2str(alpha),... 
‘) truncated to \pm’,num2str(D),’T’]) 
axis(‘equal’ ); 


Popular Symbol Alphabets 
QAM : “quadrature amplitude modulation” 


PAM : “pulse amplitude modulation” 
PSK : "phase shift keying" 


Im Im Im 


16-QAM 8-PAM 8-PSK 


Note that: "QPSK" = 4-QAM = 4-PSK and "BPSK" = 2-PAM = 2-PSK 


When the alphabet entries are spaced by A and picked with equal probability, 
the symbol variance o? = E {a [n] Pt obeys: 


alphabet M?-QAM M-PAM M-PSK 
A’ (is? A? (nf2 i 
Og" 6) (M _ 1) a (M _ 1) Asin? (1/M) 


Decision Regions 


A reasonable way to infer the transmitted symbol a|n] from the received 
sample y|n] is to decide that a|n] was the alphabet element nearest to y[n]. 


Nearest-element decision making is equivalent to using decision regions 
whose boundaries are equidistant from the two nearest alphabet elements: 


Im 


Im Im 


When a[n] = a, the symbol error rate (SER) equals the probability that y|n| 
lies outside the decision region corresponding to alphabet member a. Writing 
y|n] = a[n] + e[n], we represent the cumulative effect of noise and ISI by 
the error e|n]. Usually we model e[n] as a Gaussian random variable with 
mean 0 and variance o,”. 


Symbol Error Rate (SER) for M-PAM 


Let's first consider an M-PAM alphabet, where a[n] € R. Since the decision 
regions show that Im {y/n]} is not useful, we'll consider only the real parts 
of y[n] and e[n]. 


When a|n]| = a, we have y|n ls = a+ e[(n], implying that y/n] is Gaussian 
with mean a and variance o,”, abbreviated as “WY (a, o yr, This is illustrated 
below for the case of 4-PAM: 
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a 
Formally, we say that pyjnjjajnj (y|a), the probability density function (pdf) of 
y|n| conditioned on a[n] = . obeys 
Equation: 


Basically, Pyinjjain| (y|@) tells us how likely it is that y{n| = y given that 
aln| =a. 


Consider first the case where a is an “interior” (not an “edge”) element of the 
symbol alphabet. Given that a[n] = a, we make an error when 


y[n] < a— or when y(n] > a + 4. To find the probability of the latter 


error event, 1.e., 
Equation: 


Pr {y[n] >a+ 4 | aln] = a}, 


we integrate Pyinjjain] (y|a) over y € (a + 4, oo): 


Equation: 
: eo (y—a)” 
/ Pylnjla{n| (yla)dy = FS oP | - Ga dd 
oS at4 4/20? o 


N (a,08) 


The integral represents the shaded area below: 


\ un Py{nj\a{n] (y | a) 


\ 

7 
\ 
\l 


This integral is often solved via 
Equation: 


L ae (4 Ja o(=5*). 


— 
MN (u,07) 


using the “Q function”: 


While the Q function is not represented in MATLAB, it can be calculated 
using the “complementary error function” erfc: 
Equation: 


Q(z) = 5 erfc (=) 


In any case, the latter error event occurs with probability 
Equation: 
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By symmetry, the former error event probability is also 
Equation: 


Pr {y[n] <a—4 | a[n|=a} = o( a ). 


20 ¢ 


Since these two events are disjoint, the probability of making a decision 
error on an interior symbol equals their sum: 


Equation: 
A A A 
a(<-) +(s-) = 20(3,-). 


For edge symbols, we experience half the decision error probability, since 
there is only one decision boundary to cross. 


Finally, we average over the conditional error probabilities: 
Equation: 


Pr {error} = » Pr {error|a [n] = a}Pr {a[n] = a} 


Iphabet 
a € alphabe =1/M for alla 


(a5) ar +20 ze.) Aa 
= (ar (ze) 


Using 0? = a (M = 1), we can finally write 


12 
Equation: 
M-1 3 a? 
SER wu. = 2| —— re 
M-PAM ( uM Je Jas o 


Symbol Error Rate (SER) for M7-QAM 


With QAM, we have complex-valued y|n], a[n], e[n]. We'll assume that 
Re {e|n]} and Im {e/n]} are uncorrelated and equal variance. To calculate 
SER, we can re-use the PAM approach with a few modifications: 


1. integration is done on the complex plane, 
2. de?-variance e[n|=> % variance Re {e|n|} & Im {e[n]}, 
3. M?-QAM has 4 corner points, 4(// — 2) edge points, and 
M? — 4M + 4 interior points, 
4. calculate Pr {error|a|n] = a} via 1— Pr {correct|a|n] = a}, since 


the regions of integration are simpler: 


Pr{error|a{n] = a}: 


Pr{correct|a{n] = a bin 


After a bit of algebra, we find 
Equation: 


M-1 3 o? 


Bit Error Rate (BER) and Gray Coding 


With an M-ary alphabet, there are log, M bits per symbol, so 1 symbol error 
could cause up to log, M bit errors. 


Gray coding is a clever way of mapping bits to symbols so that neighboring 
symbols differ by only a single bit. Since the vast majority of errors occur 
when y(n] falls into a neighboring decision region, Gray coding yields BER 
~ SER. 


4-PAM: 16-QAM: 
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